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Abstract: We give a detailed account and extensions of a holographic flavor supercon- 
ductivity model which we have proposed recently. The model has an explicit field theory 
realization as strongly coupled N = 2 Super Yang-Mills theory with fundamental matter 
at finite temperature and finite isospin chemical potential. Using gauge/gravity duality, 
i. e. a probe of two flavor D7-branes in the AdS black hole background, we show that the 
system undergoes a second order phase transition with critical exponent 1/2. The new 
ground state may be interpreted as a p meson superfluid. It shows signatures known from 
superconductivity, such as an infinite dc conductivity and a gap in the frequency-dependent 
conductivity. We present a stringy picture of the condensation mechanism in terms of a 
recombination of strings. We give a detailed account of the evaluation of the non-Abelian 
Dirac-Born-Infeld action involved using two different methods. Finally we also consider 
the case of massive flavors and discuss the holographic Meissner-Ochsenfeld effect in our 
scenario. 
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Figure 1: Phase diagram for fundamental matter in thermal strongly-coupled Af — 2 SYM the- 
ory [32], with /i the isospin chemical potential, Mq the bare quark mass, M ~ 2MgA^^/^, A the 't 
Hooft coupling and T the temperature: In the blue shaded region, mesons are stable. In the white 
and green regions, the mesons melt. Here the new phase is stabilized while it was unstable in [32]. 
In this phase we find some features known from superconductivity. 



1. Introduction and summary 

Gauge/gravity duality has been used very successfully over the past years to obtain models 
with properties very similar to experimentally observed phenomena in both elementary and 
condensed matter physics. This is due to the fact that gauge/gravity duality is a powerful 
tool to compute observables of field theories at strong coupling [1]. In view of applications, 
it is particularly useful to turn on a temperature using an AdS black hole as gravity dual, as 
first proposed in [2]. In particular, gauge/gravity duality models which exhibit properties 
known from superconductivity have recently been constructed [3-22]. 

In this context we have recently proposed a model which represents a stringy realization 
of superconductivity in a relativistic framework [14]. The aim of the present paper is to give 
a detailed account of the setup and calculations involved in obtaining the results of [14], 
as well as further new results. 

Our model is based on embedding a probe of two coincident D7 branes into the AdS- 
black hole background, which extends the original correspondence to a field theory with 
flavor degrees of freedom in the fundamental representation of the gauge group, along the 
lines of [23-25]. Using this approach, the correspondence has also been applied to finite 
flavor charge densities [26-32]. A finite charge density for the flavor symmetry is obtained 
by giving a vev to the time component of the gauge field on the D7 brane probe. In [32], 
a finite SU{2) isospin density d was considered using a probe of two coincident D7 branes 
and switching on a gauge field A = A^dxaa^ on the probe. At high isospin densities 
d > dci this configuration becomes unstable against gauge field fluctuations X = a)- + ia^ 
and Y = a}- — la?, corresponding to vector mesons [33] in the field theory (see Figure 1) 
(see [34] for an early discussion of a similar instability). Note that we use capital letters 
for the background and small letters for the fluctuations. 
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The central result obtained in this model [14] is that there is a solution to the equation 
of motion for the field on the two flavor branes, which in addition to the finite isospin 
density has a further vev of the form A^dx^a^. This solution has lower free energy in the 
new phase shown in figure 1 and stabilizes this phase. This applies at least to stability 
against the fluctuations X and Y. Moreover we find that this phase has properties of a 
superconductor, such as a second order phase transition with a critical exponent of 1/2. 
The flavor current, which is analogous to the electromagnetic current, displays infinite dc 
conductivity and a gap in the frequency-dependent conductivity. Furthermore, there is 
a dual string-theoretical picture of the Cooper pairs, as well as dynamical generation of 
meson masses. 

The stabilization mechanism encountered in this context is related to the results of [35] . 
Within the Sakai-Sugimoto model, a vev of the form A^dx^a^ was considered in [36] in 
the context of p meson condensation. To obtain a p-wavc superconductor within a phe- 
nomenological AdS/QCD model, a vev of the form A^dx^a^ was first used in [7]. Note 
that in our model the equations of motion forbid a source term for A^, i. e. the breaking of 
the U{1) is spontaneous. 

The results presented here can be interpreted in two different ways. On the one 
hand, our results provide an important link between gauge/gravity with flavor and the 
recent gauge/gravity realizations of superconducting states e.g. in [3,7,13], as well as 
a realization of a superconducting current from a top-down (super-) gravity model, rather 
than in a phenomenological bottom-up AdS/QCD model. Moreover, the degrees of freedom 
in the dual field theory are identified explicitly. 

On the other hand, our holographic description can also be considered as a realization 
of a meson supcrfiuid. In the model presented here, the p mesons form a superfluid. In 
QCD like theories, p meson condensation is discussed in [37,38]. The corresponding picture 
is that the mesons move through the quark-gluon plasma without friction. Similarly, the 
absorption rate for radiation should be reduced for a specific frequency range [39]. For 
pions, similar results have been found within QCD: There are analytical results for small 
and very large isospin chemical potential using the chiral Lagrangian [40,41], as well as 
lattice results [42]. Unlike the lattice results, the results presented here can easily be 
generalized to finite baryon chemical potential. A pion superfluid was obtained within 
gauge/gravity duality in [18]. 

Let us compare our results to QCD. In QCD, the pion condensate is of course the 
natural state in isospin asymmetric matter. The condensation of a particle sets in if 
the isospin chemical potential is larger than the mass of this particle. According to this 
rule, the pions condense first in QCD since they are the Nambu-Goldstone bosons of the 
spontaneous chiral symmetry breaking and therefore the lightest particles. However the 
dual field theory which we consider in this letter is super symmetric at zero temperature 
and therefore chiral symmetry cannot be broken spontaneously. In this supersymmetric 
theory, the vector and scalar mesons have the same mass at zero temperature. Due to finite 
temperature effects, the mass of the vector and scalar mesons can become different as we 
increase the temperature. It is a priori unclear which particle will condense. In our model 
we checked that the vector mesons condense first such that the ^-meson condensation state. 
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which we consider in this paper, is the physical ground state of our system near the phase 
transition. 

In this paper we explain in detail the method used for evaluating the non-Abelian Dirac- 
Born-Infeld action involved. Moreover we give a detailed account of our results for finite 
quark mass and the derivation of the Meissner effect. Results similar to ours have also been 
obtained in [15] shortly after [14]. However there the non-abelian DBI action is considered 
only to second order in a'. Some of the features discussed here are also present in [15]: 
a second order phase transition as well as a gap in the frequency-dependent conductivity. 
However, since higher order terms in ol are not taken into account in [15], these authors 
are not able to see higher order excitations in the frequency-dependent conductivity. Here 
we identify these higher order peaks with meson resonances. 

Moreover, we give a string-theoretical interpretation of the formation of a supercon- 
ducting state of lower free energy: Essentially the condensation process corresponds to a 
recombination of strings. The detailed picture is as follows: A small isospin density cor- 
responds to strings stretching from the D7 branes to the horizon [26], the horizon strings. 
These are located close to the horizon. Increasing the isospin density, the horizon strings 
strongly charge one of the D7 branes positively near the horizon, and the other negatively. 
At a critical value, it becomes energetically favorable for the horizon strings to recombine 
into D7-D7 strings which are free to move into the bulk and therefore lower the charge 
density near the horizon. The D7-D7 strings form the superconducting condensate. Thus 
they provide a picture of strongly coupled holographic Cooper pairs. Note that for the 
construction considered the non-Abelian structure generated by the two D7-branes is es- 
sential since it provides a charged condensate. The non-Abelian forces play an important 
role for stabilizing the condensate. 

This paper is organized as follows: We begin with describing the field-theory action 
in section |2|. We then explain the dual supergravity setup in section ^ In particular, we 
contrast and compare two different approaches to evaluating the non-Abelian DBI action 
involved. These are an adapted symmetrized trace prescription and the expansion of the 
DBI action to fourth order. We frequently refer to these approaches as prescriptions. We 
then explain the string-theoretical picture in detail in section ^. Then we describe the 
thermodynamics of the phase transition in ^, and in particular show that it is second order 
with critical exponent 1/2. In ^ we evaluate the fiuctuations, again in the two different 
approaches for the non-abelian DBI action. We calculate the frequency-dependent conduc- 
tivity and show it displays the expected gap. Moreover we check that while the normal 
phase was unstable with respect to the fluctuations in X and Y, the superconducting 
phase is now stable with respect to these fluctuations. Finally we show in that our sys- 
tem displays a Meissner-Ochsenfeld effect. We conclude with a collection of remarks and 
suggestions in ^ 

2. Field Theory Interpretation 

In this paper we will focus on a 3+1 - dimensional N = 2 supersymmetric SU{Nc) Yang- 
Mills theory at temperature T, consisting of a = 4 gauge multiplet as well as Nf massive 
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N = 2 supersymmetric hypermultiplets {ipi, (pi). The hypermultiplets give rise to the flavor 
degrees transforming in the fundamental representation of the gauge group. The action is 
written down expUcitly for instance in [43]. In particular, we work in the large Nc limit 
with Nf <^ Nc at strong coupling, i.e. with A ^ 1, where A = gy^Nc is the 't Hooft 
coupling constant. In the following we will consider only two flavors, i. e. Nj = 2. The 
flavor degrees of freedom are called u and d. If the masses of the two flavor degrees are 
degenerate, the theory has a global U{2) flavor symmetry, whose overall U{1)b subgroup 
can be identified with the baryon number. 

In the following we will consider the theory at finite isospin chemical potential fi, which 
is introduced as the source of the operator 

Jq (xipa^-fotl' + (f)a^do(f) = nu-nd, (2.1) 

where n^/d is the charge density of the isospin fields, {(l)u,4'd) = and (tpuji^d) = V'- 
cr* are the Pauli matrices. A non-zero vev (Jq) introduces an isospin density as discussed 
in [32] . The isospin chemical potential /U explicitly breaks the U{2) ~ U{1)bX SU (2)/ flavor 
symmetry down to U{1)b x U (1)3, where U{1)3 is generated by the unbroken generator o"^ 
of the SU{2)i. Under the ^7(1)3 symmetry the fields with index u and d have positive and 
negative charge, respectively. 

However, the theory is unstable at large isospin chemical potential [32]. The new 
phase is sketched in figure |l|. We show in this paper (see also [14] ) , that the new phase is 
stabilized by a non-vanishing vacuum expectation value of the current component 

J3 oc ipa^'ystp + (pcr^dscj) = Tpulsipd + ipdlsi'u + bosons . (2.2) 

This current component breaks both the 50(3) rotational symmetry as well as the remain- 
ing Abelian C/(l)3 fiavor symmetry spontaneously. The rotational 5*0(3) is broken down to 
50(2)3, which is generated by rotations around the axis. Due to the non-vanishing vev 
for J3 , fiavor charged vector mesons condense and form a superfluid. Let us emphasize that 
we do not describe a color superconductor on the fleld theory side, since the condensate is 
a gauge singlet. 

Moreover, our results have a condensed matter interpretation: Our model can be con- 
sidered as a holographic p-wave superconductor. The global U{l)s in our model is the 
analog of the local C/(l)em symmetry of electromagnetic interactions. So far in all holo- 
graphic models of superconductors the breaking of a global symmetry on the fleld theory 
side is considered. In our model, the current corresponds to the electric current Jem. 
The condensate (J3) breaks the U{1)3 spontaneously. Therefore it can be viewed as the 
superconducting condensate, which is analogous to the Cooper pairs. Since the conden- 
sate (J3) transforms as a vector under spatial rotations, it is a p-wave superconductor. 
- Strictly speaking, for a superconductor interpretation it would be necessary to gauge 
the global U{1)3 symmetry which is broken spontaneously in our model. A spontaneously 
broken global symmetry corresponds to a superfluid. However, many features of super- 
conductivity do not depend on whether the U{1)3 is gauged. One exception to this is 
the Meissner-Ochsenfeld effect. To generate the currents expelling the magnetic field, the 
[7(1)3 symmetry has to be gauged. This is discussed further below in section VII. 
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3. Holographic Setup 



d^' = ^{ + /d^'l + (-) (d^?' + , (3.1) 



3.1 Background and brane configuration 

We consider asymptotically AdS^ x space-time. The AdS^ x geometry is holograph- 
ically dual to the = 4 Super Yang- Mills theory with gauge group SU{Nc). The dual 
description of a finite temperature field theory is an AdS black hole. We use the coordinates 
of [26] to write the AdS black hole background in Minkowski signature as 

with dO^ the metric of the unit 5-sphere and 

f{g) = l-^, = l + (3.2) 

where R is the AdS radius, with 

i?^ = AngsNc = 2A a'^ . (3.3) 



The temperature of the black hole given by (3.1) may be determined by demanding regu- 
larity of the Euclidean section. It is given by 

T = fL^ (3.4) 

In the following we may use the dimensionless coordinate p = q/ qh, which covers the range 
from the event horizon at p = 1 to the boundary of the AdS space at p — > oo. 

To include fundamental matter, we embed A'^^ coinciding D7-branes into the ten- 
dimensional space-time. These D7-branes host flavor gauge fields A^ with gauge group 
U{Nf). To write down the DBI action for the D7-branes, we introduce spherical coordi- 
nates {r, r^a} in the 4567-directions and polar coordinates {L,(f)} in the 89-directions [26]. 
The angle between these two spaces is denoted hy 9 {0 < 9 < tt/2). The six-dimensional 
space in the 456789-directions is given by 

d0^ + 0^dnl= dr^ + r^dnl + dL^ + L'^d4>'^ = dg^ + g^{d9^ + cos'^ 9d4>'^ + sm'^ 9dnl) , (3.5) 

where r = gsm9, = r"^ + and L = gcos9. Due to the SO {4) rotational symmetry in 
the 4567 directions, the embedding of the D7-branes only depends on the radial coordinate 
p. Defining x = cos 9, we parametrize the embedding by x = xip) and choose = using 
the S0{2) symmetry in the 89-direction. The induced metric G on the D7-brane probes is 
then 

ds\G) = -4 (-^dt' + fdxA + :^i^4±4^d,2 + R\l - x')dnl . (3.6) 



The square root of the determinant of G is given by 



G = ^g^ffii -x')\/i-x^ + g'id.xf , (3.7) 
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where /13 is the determinant of the 3-sphere metric. 

As in [32] we introduce a SU{2) isospin chemical potential /i by a non-vanishing time 
component of the non-Abelian background field on the D7-brane. The generators of the 
SU (2) gauge group are given by the Pauli matrices a*. Due to the gauge symmetry, we may 
rotate the flavor coordinates until the chemical potential lies in the third flavor direction, 



= lim Ai{p) . 

p—fOO 



(3.8) 



This non-zero gauge field breaks the SU{2) gauge symmetry down to U{1)3 generated by 
the third Pauli matrix a^. The spacetime symmetry on the boundary is still 50(3). Notice 
that the Lorentz group 5*0(3, 1) is already broken down to 50(3) by the finite temperature. 
In addition, we consider a further non-vanishing background gauge field which stabilizes 
the system for large chemical potentials. Due to the symmetry of our setup we may choose 
A^dx^a^ to be non-zero. To obtain an isotropic configuration in the field theory, this new 
gauge field ^3 only depends on p. Due to this two non-vanishing gauge fields, the field 
strength tensor on the branes has the following non-zero components. 



F. 



03 



03 



30 



7 



d,Al 



^0^3 ' 



(3.9) 



3.2 DBI action and equations of motion 

In this section we calculate the equations of motion which determine the profile of the 
D7-brane probes and of the gauge fields on these branes. An discussion of the gauge field 
profiles, which we use to give a geometrical interpretation of the stabilization of the system 
and the pairing mechanism, may be found in section ^. 

The DBI action determines the shape of the brane embeddings, i. e. the scalar fields 
(j), as well as the configuration of the gauge fields A on these branes. We consider the case 
of Nf = 2 coincident D7-branes for which the non-Abelian DBI action reads [44] 



5dbi = -Td7 Str J d^e \/ditQ 



det (Pab [E^, + E^iiQ-' - 6y^Ej,] + 2^a'F, 



ab 



(3.10) 



with 



= 6'j +i2TTa'[<^>\^'']Ekj 



(3.11) 



and Pab the pullback to the Dp-brane, where for a Dp-brane in d dimensions we have 



p, v = 0, . . . , {d — 1), a, b = 0, 
case we set p = 7, d = 10, B e 



..,p,i,j = {p+l),...,{d- 1), E^^ = g^y + B^,y. In our 
0. As in [32] we can simplify this action significantly by 
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using the spatial and gauge symmetries present in our setup. The action becomes 
Sdbi = -Tor jd^C Str^l det(G + 27ra'F)| 



-Td7 J d^e Str 

1 



where in the second hne the determinant is calculated. Due to the symmetric trace, all 
commutators between the matrices a* vanish. It is known that the symmetrized trace 
prescription in the DBI action is only valid up to fourth order in a' [45,46]. However 
the corrections to the higher order terms are suppressed by Ny'^ [47] (see also [48]). Here 
we use two different approaches to evaluate the non-Abelian DBI action ( p. 121 ). First, we 
modify the symmetrized trace prescription by omitting the commutators of the generators 
CT* and then setting (cr*)^ = 1 (see subsection 3.2.1 below). This prescription makes the 



calculation of the full DBI action feasible. Second, we expand the non-Abelian DBI action 



to fourth order in the field strength F (see subsection 3.2.2 ). 



3.2.1 Adapted symmetrized trace prescription 

Using the adapted symmetrized trace prescription defined above, the action becomes 



Sbbi = - TmNf j d^C 
TmNf 



1 + G^^G^^ (i73^)2 ^ ^33^44 ^pl^^^ ^ ^00^33 (^2^)2 

d'e £^V/(l-/)T(p,x,i), 



(3.13) 



with 



T(p,x,i) 



1 - + p\dpxf - ¥ (1 - x') (dp^if + 5(1 - x') Mr ' 



27^ 



p 



f 



(3.14) 



where the dimensionless quantities p = q/qh and A = {2Tra')A/ qh are used. To obtain 
first order equations of motion for the gauge fields which are easier to solve numerically, 
we perform a Legendre transformation. Similarly to [26, 32] we calculate the electric dis- 
placement Pq and the magnetizing field which are given by the conjugate momenta of 
the gauge fields and A^, 



pI 



6{d,Al) ■ 



(3.15) 
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In contrast to [26,27,29,32], the conjugate momenta are not constant any more but depend 
on the radial coordinate g due to the non-Abehan term ^gAg in the DBI action. For the 
dimensionless momenta ^iid defined as 

P 



P 



27ra'NfTD7Qji 



,3 ' 



we get 



PO — „ ™, T^ ' Ps — ^ 



2fT{p,x,A) 2T{p,x,A)) 
Finally, the Legendre-transformed action is given by 



(3.16) 
(3.17) 



'S'dbi = 'S'dbi - J d*^^ 
TDjNf 



J d«e Q^fm - x')v'i-x'+p'(W np^ X, A,p) , 



(3.18) 



with 

V{p,X,A,p) 



27^ 



1 + 



^6/3(1 _ ^2)3 p6//2(l_;^2)3 



(3.19) 

Then the first order equations of motion for the gauge fields and their conjugate momenta 
are 



Q^^3 ^ v^/i-x'+p\dpxy 



dpAl 



dpPl 



pIw{p,x,a,p), 



2v/l-x^ + p2(a;^)2 

p3j(l_^2)2 P3^(P' X, A, p) . 

f{l-X^)^l-x' + p'{dpX?c^ 



271'^ pfW{p,x, A, p) 



{M) Al 



(3.20) 



f(i-x')^i-x' + p'{d,xrc' (^,y^, 

2Tr^ pfWip,x,A,p) 



^3 ' 



with 



W{p,x,A,p) 



\ 



1 + 



(3.21) 



P«/-^(1-X2)« p6//2(l-x2)3 

For the embedding function x we get the second order equation of motion 



p'ff{l-X^){dpX) 
^l-x' + p2(d,xr 



V 



p^ffx 



V^-x' + pHd,xy 

24{l-x' + p\d,x?) 

/3^6(l_^2)3 



[3{l-x')+2p'{d,xf]V 



wlm'-^jApi) 



(3.22) 
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We solve the equations numerically and determine the solution by integrating the equations 
of motion from the horizon at p = 1 to the boundary at /? = oo. The initial conditions may 
be determined by the asymptotic expansion of the gravity fields near the horizon 



^0 



^3 



Pi 
Pi 



Co 



(p-i) 



+ 0{(p-lf) 



Co 



8^2 



2\3 



X = Xo 



3xo(l-Xo) 



r(p-l) 



+ 0{{p 
' + 0{ip 



1) 



I) 



I) 



(3.23) 



where the terms in the expansions are arranged according to their order in p — 1- For 
the numerical calculation we consider the terms up to sixth order in p — 1. The three 
independent parameter bo, cq and xo ™ay be determined by field theory quantities defined 
via the asymptotic expansion of the gravity fields near the boundary, 



^0 



Pi 



Pi 



X 



"0 



4 



+ 



7V4 

47r2p2 



+ o ip-') 

+ O {p~') 
+ o {p~') 

+ o ip-') 



(3.24) 



m 
P 



+- + [p 
pi 



According to the AdS/CFT dictionary, p is the isospin chemical potential. The parameters 
d are related to the vev of the flavor currents J by 



2§(J|) 



dl 



2t(ji) 



NfNcVXT^ ' NfNcVXT^ 
and m and c to the bare quark mass Mq and the quark condensate {ipil^), 



2Mn 



m 



(3.25) 



(3.26) 



y/XT ' VXNfN^T^ ' 

respectively. There are two independent physical parameters, e. g. m and p, in the grand 
canonical ensemble. From the boundary asymptotics (3.24), we also obtain that there is 
no source term for the current J3. Therefore as a non-trivial result we find that the ?7(1)3 
symmetry is always broken spontaneously. In contrast, in the related works on p-wave 
superconductors in 2 + 1 dimensions [7,8], the spontaneous breaking of the U{1)3 symmetry 
has to be put in by hand by setting the source term for the corresponding operator to zero. 
With the constraint A^\p^oo = and the two independent physical parameters, we may fix 
the three independent parameters of the near-horizon asymptotics and obtain a solution 
to the equations of motion. 
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3.2.2 Expansion of the DBI action 

We now outline the second approach which we use. Expanding the action ( |3.12| ) to fourth 
order in the field strength F yields 



5dbi = -TmNf / d«e 



2 8 



(3.27) 



where % consists of the terms with order i in F. To calculate the 7^, we use the following 
results for the symmetrized traces 



2a : Str 
4ct : Str 



0-^ 



N 



f 



Nf , Str 



3 ' 



(3.28) 



where the indices i,j are distinct. Notice that the symmetric trace of terms with unpaired 
a matrices vanish, e. g. Str[(TV-^] = NfS^K The % are given in the appendix 

To perform the Legendre transformation of the above action, we determine the con- 
jugate momenta as in ( 3.15| ). However, we cannot easily solve these equations for the 
derivative of the gauge fields since we obtain two coupled equations of third degree. Thus 
we directly calculate the equations of motion for the gauge fields on the D7-branes. The 
equations are given in ( [A.3| ) of the appendix. 

To solve these equations, we use the same strategy as in the adapted symmetrized 
trace prescription discussed above. We integrate the equations of motion from the horizon 
at p = 1 to the boundary at p = oo numerically. The initial conditions may be determined 
by the asymptotic behavior of the gravity fields near the horizon 



^0 
^3 



X 



a2(/3- 1) 



+0[{p 



I) 



XO + 



3(af + 4ai - 8)xo 
4(3a| + Aal + 8) 



(3.29) 



{p-lf +0{{p-l) 



For the numerical calculation we use the asymptotic expansion up to sixth order. As in 
the adapted symmetrized trace prescription, there are again three independent parameters 
a25&OiXo- Since we have not performed a Legendre transformation, we trade the indepen- 
dent parameter cq in the asymptotics of the conjugate momenta pg in the symmetrized trace 
prescription with the independent parameter a2 (cf. asymptotics in equation ( ^.23D ). How- 
ever, the three independent parameters may again be determined in field theory quantities 
which are defined by the asymptotics of the gravity fields near the boundary 
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(3.30) 
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Figure 2: Sketch of our string setup: The figure shows the two coincident D7 branes stretched from 
the black hole horizon to the boundary as a green and a blue plane, respectively. Strings spanned 
from the horizon of the AdS black hole to the D7-branes induce a charge at the horizon [26,27,32]. 
However, above a critical charge density, the strings charging the horizon recombine to D7-D7 
strings. These D7-D7 strings are shown in the figure. Whereas the fundamental strings stretched 
between the horizon and the D7-brane are localized near the horizon, the D7-D7 strings propagate 
into the bulk balancing the flavorelectric and gravitational, i.e. tension forces (see text). Thus 
these D7-D7 strings distribute the isospin charges along the AdS radial coordinate, leading to 
a stable configuration of reduced energy. This configuration of D7-D7 strings corresponds to a 
superconducting condensate. 



The independent parameters /x, dg, ^3, m, c are given by field theory quantities as presented 
in (3.25) and ( |3.26| ). Again we find that there is no source term for the current J3, which 
implies spontaneous symmetry breaking. Therefore the independent parameters in both 
prescriptions are the same and we can use the same strategy to solve the equations of 
motion as described below ( 3.26| ). 



4. String Theory Picture 

In this section we give a string theory interpretation, i. e. a geometrical picture, of the 
formation of a new phase, for which the field theory is discussed in section |^. We show 
that the system is stabilized by dynamically generating a non-zero vev of the current 
component J3 dual to the gauge field on the brane. Moreover, we find a geometrical 
picture of the pairing mechanism which forms the condensate (J3), the Cooper pairs. 

Let us first describe the unstable configuration in absence of the field ^3. As known 
from [26,27,32], the non-zero field is induced by fundamental strings which are stretched 
from the D7-brane to the horizon of the black hole. In the subsequent we call these strings 
'horizon strings'. Since the tension of these strings would increase as they move to the 
boundary, they are localized at the horizon, i.e. the horizon is effectively charged under 



the isospin charge given by {2A). By increasing the horizon string density, the isospin 
charge on the D7-brane at the horizon and therefore the energy of the system grows. 
In [32], the critical density was found beyond which this setup becomes unstable. In this 
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case, the strings would prefer to move towards the boundary due to the repulsive force on 
their charged endpoints generated by the flavorelectric field = Fq^ = —dgAff. 

The setup is now stabilized by the new non-zero field Al. This field is induced by D7- 
D7 strings moving in the direction. This movement of the strings may be interpreted as 
a current in direction which induces the magnetic field -63^ = = —dgA^. Moreover, 
the non-Abelian interaction between the D7-D7 strings and the horizon strings induces an 
flavorelectric field = F^q = -f/\^AlAl. 

From the profile of the gauge fields and their conjugate momenta (see figure ^) we 
obtain the following: For A^ = 0, i.e. in the normal phase (T > Tc), the isospin density 
is exclusively generated at the horizon by the horizon strings. This can also be understood 
by the profile of the conjugate momenta Pq (see figure ^ (c)). We interpret Pq{p*) as the 
isospin charge located between the horizon at p = 1 and a fictitious boundary at p = p* . In 
the normal phase, the momentum pg is constant along the radial direction p (see figure ^(c), 
blue curve), and therefore the isospin density is exclusively generated at the horizon. In 
the superconducting phase where ^3 / 0, i. e. T < Tc, the momentum Pq is not constant 
any more. Its value increases monotonically towards the boundary and asymptotes to dg 
(see figure Q(c), red curve). Thus the isospin charge is also generated in the bulk and not 
only at the horizon. This decreases the isospin charge at the horizon and stabilizes the 
system. 

Now we describe the string dynamics which distributes the isospin charge into the bulk. 
Since the field A^ induced by the D7-D7 strings is non-zero in the superconducting phase, 
these strings must be responsible for stabilizing this phase. In the normal phase, there are 
only horizon strings. In the superconducting phase, some of these strings recombine to form 
D7-D7 strings which correspond to the non-zero gauge field A^ and carry isospin charge ^. 
There are two forces acting on the D7-D7 strings, the fiavorelectric force induced by the field 
E^ and the gravitational force between the strings and the black hole. The flavorelectric 
force points to the boundary while the gravitational force points to the horizon. The 
gravitational force is determined by the change in effective string tension, which contains 
the g dependent warp factor. The position of the D7-D7 strings is determined by the 
equilibrium of these two forces. Therefore the D7-D7 strings propagate from the horizon 
into the bulk and distribute the isospin charge. 

Since the D7-D7 strings induce the fleld A^, they also generate the density ^3 dual 
to the condensate (J3), the Cooper pairs. This density ^3 is proportional to the D7-D7 
strings located in the bulk, in the same way as the density dg counts the strings which 
carry isospin charge [26]. This suggest that we can also interpret p}^{p*) as the number of 
D7-D7 strings which are located between the horizon at p = 1 and the fictitious boundary 
at p = p*. The momentum p^ is always zero at the horizon and increases monotonically in 
the bulk (see figure ^ (d)). Thus there are no D7-D7 strings at the horizon. 

The double importance of the D7-D7 strings is given by the fact that they are both 
responsible for stabilizing the superconducting phase by lowering the isospin charge density, 
as well as being the dual of the Cooper pairs since they break the U{l)s symmetry. 

^Note that the D7-D7 strings are of the same order as the horizon strings, namely Nf /Nc, since they 
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Figure 3: Profiles of the relevant dimensionless gauge fields A on the D7-branes and their dimen- 
sionless conjugate momenta p versus the dimensionless AdS radial coordinate p near the horizon 
at p = 1. The different curves correspond to the temperature T — (blue) and T O.dT^ (red). 
The plots are obtained at zero quark mass m = and by using the adapted symmetrized trace 
prescription. Similar plots may also be obtained at finite mass m ^ and by using the DBI ac- 
tion expanded to fourth order in F. These plots show the same features: (a) The gauge field Aq 
increases monotonically towards the boundary. At the boundary, its value is given by the dimen- 
sionless chemical potential p. (b) The gauge field A^ is zero for T > Tc. For T < Tc, its value 
is non-zero at the horizon and decreases monotonically towards the boundary where its value has 
to be zero, (c) The conjugate momentum pg of the gauge field A^ is constant for T > Tc- For 
T < Tc, its value increases monotonically towards the boundary. Its boundary value is given by the 
dimensionless density c?q. (d) The conjugate momentum pg of the gauge field A^ is zero for T > Tc- 
For T < Tc, its value increases monotonically towards the boundary. Its boundary value is given 
by the dimensionless density —d\. 

5. Thermodynamics & Phase Transition 

In this section we study the thermodynamics of the fundamental matter sector which is 
dual to the thermal contributions of the D7-branes. According the AdS/CFT dictionary 
the partition function Z of the boundary field theory is given in terms of the Euclidean 
on-shell supergravity action lon-sheih 

Z = e--^°-*rii . (5.1) 

Thus the thermodynamical potential, i. e. in the grand canonical ensemble the grand po- 
originate from the DBI action [27]. 



-14- 



tential J7, is proportional to the Euclidean on-shell supergravity action, 

Q = -rinZ = T/on-sheii. (5.2) 

To calculate the thermal contributions of the D7-branes Oy, we have to determine the 
Euclidean version of the DBI action ( 3.12 ) on-shell. 

5.1 Adapted symmetrized trace prescription 

Now we calculate the Euclidean on-shell action using the adapted symmetrized trace pre- 
scription defined in section First we perform a Wick rotation in the time direction. 
Next we renormalize the action by adding appropriate counterterms let (see [49] for a re- 
view). Since the action ( 3.13 ) has the same structure as the DBI action in [26,32,50] and 
since in these cases the counterterms do not depend on the finite densities, we can write 
the counterterms as in [26,32,50], 



J-Ct — T- 



( 2 
I /'max 



r 



Amc 



(5.3) 



where /3max is the UV-cutoff and 

^ TD^V^wo\{S'')NfQ\j ^ XNcMfV^T^ 
^ AT 32 ' 

with Minkowski space volume V3 . Then the renormalized Euclidean on-shell action may 
simply be written as 



(5.4) 
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G{m,ji) 



\Pmin 



m 



Amc 



(5.5) 



where pmm determines the minimal value of the coordinate p on the D7-branes, i. e. pmin = 1 
for black hole embeddings which we consider exclusively in this paper, and 



G{m,fl) 



dpip'ff{l-x^)T{p,x,A)-{p' 



pm) 



In figure |^ we plot the dimensionless grand potential W7 defined as 
We also calculate the specific heat Cj 



TIr = i-^^ Wj . 

^ 32 ' 



C7 



(5.6) 



(5.7) 



(5.8) 



where Fj is the free energy in the canonical ensemble ^. Using (3.26), the expression for 
the specific heat can be rewritten as 



C7 



^ I VlT-j - 6m- ^ 



32 



dm 



m' 



dmi? 1 



(5.9) 



^Note that in [14] we considered the second derivative of the grand potential in the grand canonical 
ensemble. Although not simply related to the specific heat, that quantity also shows a discontinuity 
indicating the second order phase transition. 



-15- 



with F7 = {XN f NcVsT'^ / 32) T-j. In figure ^ we plot the dimensionless specific heat C7 
defined as 



XNfNrT^ 



7 



32 



C7 



(5.10) 



Below we discuss that our results hint towards a boson liquid being present. We also 
comment on the relation of our results to [51]. 

5.2 Expansion of the DBI action 

The thermodynamical behavior in the case of the expanded DBI action may be determined 
in the same way as for the adapted symmetrized trace prescription. As mentioned above, 
the counterterms are needed to regularize the action and do not depend on the gauge 
fields on the brane. They are needed to regularize the divergence which occurs due to 
the infinite volume of the AdS space. For vanishing gauge fields on the brane, the action 
calculated using the adapted symmetrized trace prescription and the action expanded in 
field strength F coincide. Thus the counterterms discussed above may also be used to 
regularize the action expanded in F. For the renormalized Euclidean action we obtain 

,2n2 



^ = G(m,/i)-i 
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vrmin 



m 



r 



where 



G{m,fl) 



dp 



-G 



2 8 



{p^ — prn) 



(5.11) 



(5.12) 



with T2 and 74 as defined in equations (|3.271 ) and (A.l). In figure g we also plot the 
dimensionless grand potential W7 obtained by the action expanded to fourth order in F 
and defined as 

XNfN,V^T^ 



TIr 



32 



-W7. 



(5.13) 



Note that we have not computed the Legendre-transformed action for the DBI action 



expanded to fourth order due to technical complications discussed in section |3.2.2| . Thus 
we do not have direct access to the quantities in the canonical ensemble, such as the free 
energy needed to compute the specific heat. A calculation exploiting thermodynamical 
relations between the ensembles is postponed to future work. 

5.3 Thermodynamical consistency 

In the following we confirm that both definitions of the grand potential, ( ^.7] ) and ( 5.13| ), 
give consistent thermodynamics. For this purpose, we have to vary the grand potential with 
respect to the thermodynamical variables. In the gravity setup this variation is induced by 
a variation of the gravity fields. Here we only consider a variation of the gauge field A^, 
since the variation of the embedding function x is the same as in [26]. Using the equations 
of motion, the variation of the dimensionless grand potential with respect to the gauge 
field A^j^ reduces to a boundary term. 
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(5.14) 
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Figure 4: The dimensionless grand canonical potential calculated using the adapted sym- 
metrized trace prescription versus temperature at zero quark mass Mq = 0: Below T = Tc the 
superconducting phase (red line) is thermodynamically preferred over the normal phase (blue line). 



Here C is given by the integrand of the grand potential. Using the asymptotic behavior 
of the gauge fields ( ^331) , (|3^ ), ( ^ ) (|]2|), the contribution of the fields Al and p\ 



vanishes and the result coincides with [26], 

SW7 = -pI{oo)5AI{oo) = -dlSfi . (5.15) 

Therefore we confirm that is the thermodynamical potential in the grand canonical 
ensemble, namely the grand potential. We also see that the density ^3 is not a thermody- 
namical variable. 

5.4 Results & Comparison of the two prescriptions 

We now compare the results obtained by evaluating the non-Abelian DBI action using the 
adapted symmetrized trace prescription or the expansion to fourth order. 

Figure ^ shows that for non-zero gauge field ^3, a phase transition occurs. The pres- 
ence of this phase transition does not depend on the prescription used, as seen from the 
comparison of the potential derived using the adapted symmetrized trace prescription (see 
figure |5|) with the results derived from the expanded DBI action (see figure |^ . The transi- 
tion is second order since there is a discontinuous step in the specific heat C7, see figure |9|. 

The temperature scale in the massless case is defined by 



(5.16) 







where jlc ~ 2.85 and ((io)c ^ 20.7 as obtained from the adapted symmetrized trace pre- 
scription (see section |3. 2. ID . For the expanded DBI action we equivalently obtain flc ~ 2.48 
and {dl)c « 9.03. In the massive case we fix a quark mass Mq, then a chemical potential 
fj, measured in units of Mg. Thus the temperature scale is defined as 

T rrir 

7^ = — , 5.17 

Tc m 
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Figure 5: Same as figure || close to the phase transition. 
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Figure 6: Grand canonical potential computed from the expanded DBI action at vanishing quark 
mass Mq = 0. The qualitative behavior is in agreement with the results from the adapted sym- 
metrized trace prescription (see figure ||) . 



where rric is the value of the mass parameter at which the phase transition occurs. 

From the behavior of the order parameter dg, see in figure ^, near the critical temper- 
ature, we obtain the critical exponent of the transition to be 1/2, which coincides with the 
result obtained from Landau theory. It is also determined numerically from a fit within 
10 percent accuracy, see blue dashed curve in figure 0. This exponent is independent from 
the quark mass as seen from the black curve at fJ^/Mq = 3 shown in figure 0. We have 
also confirmed numerically that the critical exponent of the order parameter does not 
depend on the prescription (adapted symmetrized trace or expanded DBI). Since the order 
parameter and the density increase rapidly as T — > 0, we expect that the probe 
approximation breaks down near T = 0. 

We define the density of superconducting charges dg in the following way 



pUpb) - pUph) ^ dp -Co 
pUpb) dl 

This identification was suggested in [7]. Beyond the arguments given there, we have the 



4 = ""^""4, ""^""^ =^^Sr^- (5.18) 

Po(PB) dg 
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Figure 7: The order parameter d\ defined in ( 3.2'j ) and obtained from the adapted symmetrized 
trace prescription versus temperature T: The case of vanishing quark mass (red curve) shows the 
same behavior near as that at finite mass where n/Mq = 3 is fixed. In both cases the order 

1/2 



parameter vanishes with a critical exponent of 1/2. This is visualized by the fit 55(1 
(dashed blue curve). 



Figure 8: Superconducting density ds = (dp — co)/do versus temperature T: In both, the massless 
(red curve) and the massive case at ^i/Mq — 3 (blue curve), the superconducting density ds vanishes 
linearly at the critical temperature. This is visualized by the fit 6.8(1 — T/Tc) (dashed blue curve). 



advantage of exactly knowing the dual field theory interpretation and having a string 
picture, see section |^. The numerator of equation ( ^.18| ) counts the isospin charges present 
in the bulk excluding those localized at the black hole horizon Pq{ph) = cq- Following the 
ideas described in section ^ these charges may be identified with the D7-D7 strings in the 
bulk, which are dual to the Cooper pairs. As expected from Ginzburg-Landau theory, this 
superconducting density dg vanishes linearly at Tc in massive and massless cases regardless 
of the prescription chosen, see figure ^. 

In figure |9| the dimensionless specific heat computed from the adapted symmetrized 
trace prescription for massless quarks is shown versus the temperature. It clearly shows 
a discontinuity at the phase transition, stressing that the observed transition is second 
order. Note that the dimensionful specific heat is given by C-j oc T'^Cj, where C7 is the 
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Figure 9: The flavor brane contribution to the specific heat as computed from the adapted sym- 
metrized trace prescription in the massless case. The blue hne corresponds to the normal phase 
with A\ = 0, while the red one corresponds to the superconducting phase with A3 ^ 0. Note that 
the total specific heat is always positive although the flavor brane contribution is negative. The 
divergences near T = in both phases can be attributed to the missing backreaction in our setup. 
We read of from the numerical result that near the critical temperature, the dimensionless specific 
heat is constant in the superconducting phase. This implies that the dimensionful specific heat is 
proportional to T^. This temperature dependence is characteristic for Bose liquids. 

dimensionless specific heat shown in figure ^. 

The blue curve in figure ^ corresponds to the normal phase while the red one shows the 
behavior of the superconducting phase. A general feature of superfluids, including charged 
superfluids as in the present case, is the behavior of the specific heat near T = 0. For Bose 
liquids the heat capacity has cubic temperature dependence. In Fermi liquids, a linear 
temperature dependence dominates near T = 0. In our numerical results (figure ^ we 
see a dramatic change of specific heat in the normal phase (blue curve) in sharp contrast 
to the constant behavior, corresponding to C7 oc T^, of the superconducting phase (red 
curve) close to the critical temperature. The divergent behavior close to T = is due to 
the probe approximation which we use. Although our numerics do not allow to extend the 
calculation up to T = 0, it is tempting to extrapolate the observed behavior towards T = 
and to conclude that our setup in the superconducting phase is a Bose liquid. 

Recently a new type of quantum liquid has been discussed in [51], where the authors 
consider a brane setup at zero mass in the limit of infinite baryon density, i. e. at zero 
temperature. In that case the heat capacity was found to be proportional to T^, neither 
typical for a Fermi nor for a Bose liquid. 

In the following we argue that our result found above suggests to consider the possibility 
of a phase transition also in the baryonic setup [51]. The baryonic setup of [51] may 
be related to our normal phase by a symmetry between pure baryonic and pure isospin 
configurations found in [32]. Thus the specific heat in the baryonic setup and in our 
normal-conducting phase coincide. Since in our isospin case there is a phase transition 
to the superconducting phase which shows the usual Bose liquid behavior, we suggest the 
possibility of a similar phase transition in the baryonic setup of [51]. The possible new 
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phase may lead to the usual behavior of a Bose or Fermi liquid. 




Figure 10: Phase diagram for fundamental matter with mass m — 2Mq/{^/XT) obtained with 
the adapted symmetrized trace prescription: The blue, white and green regions are the same as 
in hgure 0, but with the unstable normal phase replaced with the superconducting phase. The 
dotted curves correspond to lines at finite mass, i.e. constant xo- These curves are parametrized 
by the density c?q. Along the blue curves the field is zero while along the red ones the field 
is non-zero. The endpoints of the red curves determine the second order phase transition to the 
superconducting phase. The dotted, red curves diverge inside the superconducting phase since the 
backreaction of the condensate on the background is not considered. This divergence determines 
the boundary of the orange region which is not reachable without backreaction. 



Summarizing our thermodynamical results in a phase diagram, we obtain figure 10 
The choice of the calculational method - adapted symmetrized trace or expansion of the 
DBI to fourth order - does not change the qualitative structure of the phase diagram. 
The blue phase indicates the known region of stable mesons surviving the deconfinement 
transition. It is separated from the white meson melting region by the meson melting 
transition (blue line), see [25-27,29,32,50,52,53]. Above a critical isospin density marked 
by the green line, a flavor-superconducting phase forms. At even higher isospin density, our 
approach, which does not include the back-reaction of the D7-brane, gives diverging order 
parameters, signalling the breakdown of this approach. This particular region is indicated 
by orange color. Note that this behavior is well-known to occur in systems without back- 
reaction. In our case, the gauge fields on the D7-brane grow arbitrarily large close to 
zero temperature. This behavior will be cured by including the back-reaction, i. e. the 
contribution of the gauge fields to the total energy- momentum tensor. For the Abelian 
Higgs model, it is shown in [13] that similar divergences are removed when these terms are 
included. 



6. Fluctuations 

The full gauge field A on the branes consists of the field A and fiuctuations a, 

i = Ajjr^dt + A^rMxs + ^r-^dx^ , (6.1) 
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where r° are the SU{2) generators. The hnearized equations of motion for the fluctuations 
a are obtained by expanding the DBI action in a to second order. We wih analyze the 
fluctuations and X = 02 + ia^) Y = — ia2- 

Including these fluctuations, the DBI action reads 



S = -Tt j d^C Str^det [G + {2^a')F] , 
with the non-Abelian field strength tensor 



rpa 



pa I 77a 



where the background is collected in 



and all terms containing fluctuations in the gauge field are summed in 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



Index anti-symmetrization is always defined with a factor of two in the following way 
\K\ = (df^A, - d,A^)/2. 

6.1 Adapted symmetrized trace prescription 

In this section we use the adapted symmetrized trace prescription to determine the fluc- 



tuations about the background we discussed in section 3.2.1. To obtain the linearized 
equations of motion for the fluctuations a, we expand the action ( |6.2| ) to second order in 
fluctuations, 



5(2) 



T7 / d^e Str 



— (2ira') I — -7^,,,, „ 



(27raO 



/\2 



-QG''^' F^,,G'"'' F, 



+ 



{2-71 a' 



,/^2 



(6.6) 



As in the previous paper [33] , we collect the metric and gauge field background in the tensor 
Q = G + {27ra')F. Using the Euler-Lagrange equation, we get the linearized equation of 
motion for fluctuations ai in the form 



=aAStr 



+ A)i + (27rQ')(g'^['^g^]^ + ^a'^^at"^!)^.;.^} 



Str 



abda 



(6.7) 



Note that the linearized version of the fluctuation field strength used in equation ( |6.7| ) is 
given by 

= 25[^<j + ^r^^iA'^al + a\Al) + O^a") . (6.8) 
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In our specific case the background tensor in its covariant form is given by 

T /|3/il. 



Q^u = G^.t'' + {2^a') (2d,Al5^^^6,^^T^ + 2d,A%y^5,^^T^ + 2J^AIA\5oi^5,^:,t'') . (6.9) 

Inversion yields the contravariant form needed to compute the exphcit equations of motion. 
The inverse of Q is defined as Q^^'^Q^^i = 5^,t^ ^. The non-zero components of Q^'^ may be 



found in the appendix B.l 



Fluctuations in 02 : For the fluctuation with zero spatial momentum, we obtain the 
equation of motion 



(a^r + ^(ai)' 



with 
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4 = 0, 



(6.10) 
(6.11) 



Fluctuations in X = a2 + la^, Y = a2 — ia^: For the fluctuations X and Y with zero 
spatial momentum, we obtain the coupled equations of motion 
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where the component of the inverse background tensor may be found in the appendix B.l 
index symmetrization is defined Q^^^^ = [Q^^ + Q^^)/2 and 



3- ^ i3. 



2^/2n 



(6.14) 



^We calculate the inverse of G by ignoring the commutation relation of the t's because of symmetrized 
trace. It is important that t°t'' must not be simplified to e'^^'^r'^ since symmetrisation is not the same. 
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6.2 Expansion of the DBI action 

In this section we determine the equation of motion for the fluctuation in the background 
determined by the DBI action expanded to fourth order in F (see section p. 2. 21 ). To obtain 
the quadratic action in the field a\, we first have to expand the DBI action ( |6.2[) to fourth 
order in the full gauge field strength F, and expand the result to second order in Og. Due 
to the symmetries of our setup, the equation of motion for the fiuctuation at zero spatial 
momentum decouples from the other equations of motion, such that we can write down 
an effective Lagrangian for the fiuctuation a|. This effective Lagrangian is given in the 



appendix B.2. The equation of motion for 02 with zero spatial momentum determined by 
the Euler-Lagrange equation is given by 
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where 



n 



(6.16) 



We introduce the factors Ti^^ which may be found in the appendix |B.2| to emphasize the 
similarity to the equation of motion obtained by the adapted symmetrized trace prescription 

derol). 



6.3 Result & Comparison of the two prescriptions 
6.3.1 Conductivity 

We calculate the frequency-dependent conductivity cr{uj) using the Kubo formula. 



a(uj) = -G"'{LO,q = 0) 



(6.17) 



where is the retarded Green function of the current dual to the fiuctuation a^, which 
we calculate using the method obtained in [54] . The current Jf is the analog to the electric 
current since it is charged under the f/(l)3 symmetry. In real space it is transverse to the 
condensate. Since this fluctuation is the only one which transforms as a vector under the 
50 (2) rotational symmetry, it decouples from the other fluctuations of the system. 

The real part of the frequency-dependent conductivity Recr(u;) is presented in flgure 11 
12, 13, |l^. It shows the appearance and growth of a gap as we increase the condensate ^3 



The conductivity gap originates in a pseudo gap already present right above Tc, as can be 



seen for example from the red curve in flgure 11. By pseudo gap we mean a well-deflned gap 
in the conductivity at low frequency in which the conductivity is not identically zero [8]. 
Both calculational prescriptions - the adapted symmetrized trace prescription and the 
expansion of the DBI action to fourth order in the fleld strength - yield qualitatively very 



similar results as can be seen by comparing flgures 13 and 15. As a further distinct effect 



when using the adapted symmetrized trace prescription, flgure 11, 12, 13 show prominent 
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tt) 

Figure 11: Real part of conductivity, Re cr, in units of NfNcT/{16TT) versus the dimensionless 
frequency in — u}/{2ttT) for massless quarks computed from the adapted symmetrized trace pre- 
scription. Distinct curves correspond to T/Tc = oo (black), l(red), 0.5 (orange) and 0.28 (brown). 
By decreasing the temperature below the critical one, a gap where the conductivity is approximately 
zero appears which is a characteristic feature of a superconductor. In addition prominent peaks 



peaks which we interpret as mesonic excitations below. Increasing the quark mass Mq 
from zero to a finite value, these meson peaks become sharper, i.e. more quasiparticle-like. 
This is reminiscent of results for condensed matter systems where prominent quasiparticle 
peaks appear (e.g. [55]). In figure ^ at fixed n/Mq = 3, the sharp resonances are also 
present inside the gap. In contrast to the adapted symmetrized trace prescription, from 



the expanded DBI action we obtain less prominent peaks as seen from figure 15. In the 
conductivity obtained from the DBI action expanded to fourth order, the peaks do not 
appear until we approach small temperatures. We expect that the terms higher order in 
the field strength dominate the generation of the peaks and therefore the generation of the 
meson mass. If the higher order terms discussed here are included, we presume that the 
quasinormal modes which generate these peaks move closer to the real axis. 

Using the Kramers-Kronig relation, which connects the real and imaginary part of the 
complex conductivity, we find a delta peak at a; = in the real part of the conductivity. 
Re o"(u;) ~ ims5{u)). The corresponding n^/w-behavior in the imaginary part is visualized 
in figure |l^ As expected from Ginzburg-Landau theory, our numerics show that the 
superconducting density vanishes linearly at the critical temperature, Ug oc (1 — T/Tc) for 
T TuTc- This field theory definition of the superconducting density yields a quantity with 



the same linear scaling near Tc as found in our bulk definition ( 5.18| ) of the superconducting 



density ds- This confirms that these two quantities may be identified. Our numerics indeed 
indicate that dg and Ug are proportional to each other. 

Note that for translation invariant systems at finite density, there is a delta peak in the 
real part of the conductivity even in a normal conducting phase since the charge carriers 
cannot lose their momentum. This peak is called Drude peak. In our system, however, the 
charge carriers can dissipate their momentum although our system is translation invariant 
[56,57]. The adjoint degrees of freedom can transfer momentum at order N'^ while the 
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Figure 12: Real part of conductivity, Re cr, in units of NfNcT/{16TT) versus the dimensionless 
frequency tr — u;/(2TrT) for massive quarks at fixed fJ-/Mq = 3 computed from the adapted sym- 
metrized trace prescription. Distinct curves correspond to T/Tc = 10 (black), l(red), 0.6 (green) 
and 0.5 (orange). As in figure |ll| a gap appears. Moreover, we observe a prominent peak inside the 
gap where the conductivity is approximately zero in the green and orange curves. 
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Figure 13: Scaled real part of conductivity Re cr/(47rrD) in units of NfNcT/{16Tr) versus the 
dimensionless frequency tn = u}/{2'kT) for massive quarks at fixed fi/Mq = 3 computed from 
the adapted symmetrized trace prescription. Distinct curves correspond to T/Tc ~ 10 (black), 
l(red), 0.5 (orange) and 0.28 (brown). This figure has been scaled to asymptote to a constant in 
order to show similarity to the lower dimensional cases computed from AdS^ and to show common 
asymptotics. 

fundamental degrees of freedom only at order Nc- The adjoint degrees of freedom effectively 
act as a heat sink into which the flavor fields can dissipate their momentum. Thus we do 
not observe a Drude peak in our system. 

6.3.2 Spectral functions 

In order to identify the prominent peaks in the conductivity in terms of field theory quan- 
tities as meson resonances, we examine the spectral functions TZ = — 2ImG^ which are 
related to these conductivities by ( |6.17| ). Due to this relation the prominent peaks in the 
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Figure 14: Scaled imaginary part of conductivity ir Im cr in units of NfNcT/{16TT) versus the 
dimensionless frequency tn = uj/{2itT) for massless quarks computed from the adapted symmetrized 
trace prescription. Distinct curves correspond to T/Tc = oo (black), 1 (red), 0.9 (blue), 0.6 (green), 
0.5 (orange) and 0.28 (brown). This figure has been scaled to asymptote to a constant at Id = 0. 
This constant determines the superconducting density n^. 
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Figure 15: Scaled real part of conductivity Re(T/(47rrD) in units of NfNcT/{16Tr) versus the 
dimensionless frequency to = uj/{2ttT) for massless quarks computed from the expanded DBI action. 
Distinct curves correspond to T/Tc = oo (black), l(red), 0.6 (green), 0.5 (orange) and 0.39 (purple). 
This figure has been scaled to asymptote to a constant in order to show similarity to the lower 
dimensional cases computed from AdSi and to show common asymptotics. 



conductivity correspond to the prominent peaks in the spectral function, see figure 16. As 
known from [33, 58] the resonances appearing in the spectral functions of vector fields in 
the bulk correspond to vector meson excitations in the dual field theory. Moving to higher 
quark mass parameter m = 2.842 and fl = 3.483 near the meson melting transition, we 
compare the resulting spectrum, see figure |l^ to the supersymmetric mass formula ob- 
tained in [24]. The prominent peaks clearly approach the supersymmetric line spectrum 
from above. The same behavior was found in figure 7 of [33] , where vector mesons were con- 
sidered also close to the meson melting transition. Note that the accuracy of our numerics 



-27- 



becomes insufficient at larger frequencies around tn ~ 15 




Figure 16: Tlie spectral function TZ in units of NfNcT'^/8 versus the dimensionless frequency 
ro = uj/{2ttT) for massless quarks Mq = 0. The distinct curves corresponds to different temperature 
(for color coding see figure When decreasing the temperature, peaks which correspond to 

mesonic quasiparticles emerge. 



6.3.3 Comments on stability 

As shown in figure |l^, our setup is stable with respect to the fluctuations X and Y. 



Furthermore, figure 18 shows that the quasinormal modes of higher excitations n > 1 move 
to larger frequencies and closer to the real axis. This corresponds to the formation of stable 
massive mesons. Such a behavior is known in gauge/gravity duality for mesons which are 



*The frequency range of validity of our fluctuation numerics increases with quark mass. However, the 
numerics for the background become worse at large quark masses. 




Figure 17: Finite temperature part of the spectral function TZ — TZq with TZq = 47rrD^in units of 
NfNcT'^/8 versus the dimensionless frequency fo — llj/{2ttT) at finite mass m = 2.842 and chemical 
potential fl = 3.483. The grey lines correspond to the supersymmetric mass spectrum calculated 
in [24]. 
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built from massive quarks (e.g. [33]). Thus we observe a dynamical mass generation for 
the mesons as described in the following subsection. 

Im ro 







n = 1 ) / 


^ \^ ( n = 1 


/ / n = 


n = 0\\ 




n = \ 



Figure 18: Movement of quasinormal modes under changes of the temperature T: The different 
colors indicate the different fluctuations X (red), Y (green) and a\ (blue). The higher excitations 
of the fields X and Y behave as excitations with a non-zero quark mass. This indicates a dynamical 
generation of the meson mass. 



6.3.4 Dynamical mass generation 

In this section we discuss dynamical mas generation on the field theory and on the gravity 
side. For this issue, it is important to distinguish between the cases where the broken 
symmetry is global or local in the boundary field theory. 



Field theory observation As explained in section |3.2| , the superconducting condensate 
breaks the U{\)-i symmetry spontaneously in the field theory living on the AdS boundary. 
According to the Goldstone theorem this generates one massless Nambu-Goldstone boson 
in the boundary field theory. Let us now discuss the two different cases: 

Local C/(l)3 symmetry: If the broken symmetry is gauged, the Nambu-Goldstone boson 
is eaten by the gauge field charged under the spontaneously broken ?7(f)3 sym- 
metry. In conventional superconductors this mechanism gives mass to the photons 
which implies the Meissner-Ochsenfeld effect. 

Global ?7(I)3 symmetry: If the broken symmetry is global, there is no dynamical gauge 
field which eats the Nambu-Goldstone boson. Thus the Nambu-Goldstone boson 
remains present in the spectrum. Since this is the case in our setup, we need to iden- 
tify this Nambu-Goldstone boson in the spectrum. In general, the Nambu-Goldstone 
boson corresponds to the phase of the condensate which parametrizes the the coset 
space U{l)z/'^2- In our setup, the fluctuations X and Y defined above ( 6.12| ) are 



charged under the U{1)^ symmetry such that the Nambu-Goldstone boson can be 



found in these fluctuations (see fig. 18). 



Although in our case the broken symmetry is global, we nevertheless observe dynam- 
ical mass generation as see for instance in figures ^ and 18. Therefore in our setup, a 



more subtle mechanism then the ordinary Higgs mechanism generates the meson masses 
dynamically. Let us explain this mechanism in the dual gravity setup. 
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Gravity explanation According to the AdS/CFT dictionary, a global symmetry in the 
boundary field theory is mapped to a local symmetry in the gravity theory. Thus the 
SU{2) symmetry of the two D7-branes is local in the gravity theory. In this paper we find 
new solutions to the equations of motion which non-vanishing gauge field components 
and ^3. Since the gravity action is unchanged and thus still SU (2) invariant, this non-zero 
components break the SU{2) gauge symmetry spontaneously which generates in total three 
Nambu-Goldstone bosons. They are eaten by the gauge fields which then become massive. 
Let us now discuss the corresponding mass terms in more detail. 

First we consider a non-zero SU{2) isospin chemical potential Aq which generates 
the mass term {Alf{a^'^f for the fluctuations a^'^. This new term shifts the quadratic 
fluctuation term in the action, i.e. its energy by the value of the field Aq. For the spectral 
function of the dual field theory this results in a shift of the meson resonance peaks, see [33]. 

Second we consider the non-zero gauge field component A3 which generates a mass term 
(^13)^(0^)^ for the fluctuation a^. We expect that this mass term shifts the quasinormal 
frequencies of this system in the complex plane. That results in generation and shift of the 
resonance peaks in the spectral function as well as generation of the gap structure in the 
conductivity. To make more precise statements, the spectrum of quasinormal modes needs 
to be studied in presence of these flelds. 



7. Meissner-Ochsenfeld effect 



The Meissner effect is a distinct signature of conventional and unconventional superconduc- 
tors. It is the phenomenon of expulsion of external magnetic fields. An induced current in 
the superconductor generates a magnetic field counter-acting the external magnetic field H. 
In AdS / CFT we are not able to observe the generation of counter-fields since the symme- 
tries on the boundary are always global. Nevertheless, we can study their cause, i.e. the 
current induced in the superconductor. As usual [4-6, 9, 13, 18] the philosophy here is to 
weakly gauge the boundary theory afterwards. 

In order to investigate how an external magnetic field infiuences our p-wave supercon- 
ductor, we have two choices. Either we introduce the field along the spatial z-direction 
H^T^ or equivalently one along the y-direction i.e. -fff''"^- Both are "aligned" with the 
spontaneously broken C7(l)3-flavor direction. 

As an example here we choose a non- vanishing H^t^- This requires inclusion of some 
more non- vanishing fleld strength components in addition to those given in equation (^.S|). 
In particular we choose AQ{g,x), A^{q,x) and A2{x) = yielding the additional com- 
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F 
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2 

23 



-F' 



-Fl = diAl , 

. ^ ,/|3/ll 
l3 



(7.1) 



10 — 01 — ^1^0 ■ 

Recall that the radial AdS-direction is identified by the indices either g or 4. Amending 
the DBI-action ( p.l2| ) with the additional components (7.1), we compute the determinant 
in analogy to equation (|3.12 ). We then choose to expand the new action to second order in 
F, i.e. we only consider terms being at most quadratic in the fields. This procedure gives 
the truncated DBI action 



5dbi = - TdjNj j d«C 



1 + l^^f g00g33(f23)2 + G^^G^\Fl^f + G0°G44(f3 )2 



+ 



[G^'^fiFlf + iG'^f{Fi,f + G^'^G'HFfof + (^33)2(^3^)2) 



+ 



(7.2) 



Respecting the symmetries and variable dependencies in our specific system, this can be 
written as 



'S'dbi 



TojNf / d»^ V-G 



1 + lf^G'"'G'''\diA'oy + G''''G''\d,A'^)'' + G'"'G''\d,A'oj 



+ (G33)2(ajiJ)2 + {G'^)\Hl)^ + g00g33 g^V .^pi^2 



(27ra')2A 



with the convenient redefinitions 

A = A, X = qhx , rl- 

QH 

Rescaling the x-coordinate once more 



2TTa'Hl , 



QH P- 



V 27ra'\/A ' 

the equations of motion derived from the action (|7.3D take a simple form 



(7.3) 
(7.4) 

(7.5) 
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2^/27r G44 
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2^2^ G44 
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(7.6) 



(A3)2ii 



27r2 G44 



3 • 



^ Close to the phase transition, it is consistent to drop the dependence of the field //f on q. Away form 
the phase transition the q dependence must be included. From the boundary asymptotics it will be possible 
to extract the magnetic field and the magnetization of the superconductor. 
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Here all metric components are understood to be evaluated at R = 1 and g ^ p- 

We aim at decoupling and solving the system of partial differential equations (^^) by 
the product ansatz 

Al{p,3:)=g{p)f{x). (7.7) 

For this ansatz to work, we need to make two assumptions: First we assume that Aq is 
constant in x. Second we assume that ^3 is small, which clearly is the case near the 
transition T ^ Tc- Our second assumption prevents from receiving a dependence on x 
through its coupling to A\{p^x). These assumptions allow to write the second equation in 
(O) as 



7^1 dl^{x)-x^m 



(7.8) 



2^2^ f(x) 



All terms but the last one are independent of x, so the product ansatz ( [T.T] ) is consistent 
only if 

MMf^ = C, (7.9) 
f(x) 

where C is a constant. The differential equation (^^) has a particular solution if C = 
— (2n + 1) ,n G N. The solutions for f(x) are Hermite functions 

which have Gaussian decay at large \x\ ^ 1. Choosing the lowest solution with n = and 
Hq = 1, which has no nodes, is most likely to give the configuration with lowest energy 
content. So the system we need to solve is finally given by 



= ^P^o + ""'ISooG^r^ d,Al^^^AlMi),ip)f , 
- a g H —dpQ - 77-0 7^4 (^^0) - o /?r_^44 • 



1^(^33(^44 27r2G44^ 'J'*' 2^/27r 

Asymptotically near the horizon the fields take the form 

Al = +a2{p-lf +0{{p-l 



(7.11) 



h (7.12) 

=bo +'-^{p-lf +0{{p-ir'), 



while at the boundary we obtain 



P^ 

d. 



Al =p +^ +0{p-') 



1 

3 1 /n/^-4\ 



(7.13) 



P^ 
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T 

Figure 19: The line of critical magnetic fiefd versus critical temperature. Below this line the 
external magnetic field coexists with the superconducting condensate. Above the line the super- 
conducting condensate vanishes. We set the spatial position arbitrarily to i = 0.1 since the critical 
line does not depend on x. 



We succeed in finding numerical solutions g(/j) and A^{p) to the set of equations ( 7.11| ) 
obeying the asymptotics given by equations ([7.12D and ( 7.13| ). These numerical solutions 
are used to approach the phase transition from the superconducting phase by increasing the 
magnetic field. We map out the line of critical temperature-magnetic field pairs in figure 



19. In this way we obtain a phase diagram displaying the Meissner effect. The critical line 
in figure [l9| separates the phase with and without superconducting condensate ^3. 

We emphasize that this is a background calculation involving no fluctuations. Com- 
plementary to the procedure described above we also confirmed the phase diagram using 
the instability of the normal phase against fluctuations. Starting at large magnetic field 
and vanishing condensate ^3, we determine for a given magnetic field the tempera- 
ture Tc(H^) at which the fluctuation af becomes unstable. That instability signals the 
condensation process into the superconducting phase. 

The presence of the coexistence phase below the critical line, where the system is still 
superconducting despite the presence of an external magnetic field, is the signal of the 
Meissner effect in the case of a global symmetry considered here. If we now weakly gauged 
the flavor symmetry at the boundary, the superconducting current J3 would generate a 
magnetic field opposite to the external field. Thus the phase observed is a necessary 
condition in the case of a global symmetry for finding the standard Meissner effect when 
gauging the symmetry. 



8. Discussion and Outlook 

We found a holographic realization of superconductivity in the context of gauge/gravity 
duality with flavor at finite isospin chemical potential, for which the field theory action 
is known explicitly. The condensation process corresponds to a recombination of strings 
which leads to a thermodynamically favored configuration. On the field theory side, this 
new state may be interpreted as a p meson superfluid. We show that the superconducting 
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phase persists at finite quark mass and find a holograptiic Meissner-Odisenfeld effect. 
From a technical point of view, we compare two different approaches for evaluating the non- 
Abelian DBI action involved (expansion to fourth order and a new modified trace evaluation 
procedure which we propose) and find the same physical results in both approaches, at least 
qualitatively. 

For the future, we expect that a holographic s-wave superconductor may be constructed 
by applying the methods presented to the transverse scalars in similar holographic setups. 
It will also be interesting to apply our methods to the configuration discussed in [48] in 
order to find a purely fermionic superconductor. 
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A. Expansion of the DBI action to fourth order 



In this section we give some lengthy equations needed to calculate the expanded DBI action 
and the resulting equations of motion (see section 3.2.2| ). 

The % which determine the terms in expanded DBI action a.t i — th order in the field 
strength are given by 
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Using the Euler-Lagrange equation the resulting equations of motion are 
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B. Fluctuations 



B.l Adapted symmetrized trace prescription 



In this section we present the elements of the inverse background tensor Q^'^ defined in 



section 3.1. The diagonal elements are given by 
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l + G33g44(g^^l)S 

N 
iV 



^33 _ ^33 1 + G°°G^^((9p^o)2 



(B.l) 



-35- 



where the other diagonal elements are just Q^^ = G**. The off-diagonal elements read 
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The denominator N is given by 



N = l + G'^G^^{d,Alf + G^'G'^{d,Alf - -^{AlAlf 



(B.3) 



B.2 Expansion of the DBI action to fourth order 



In this section we calculate the quadratic action for the fluctuations using the expanded 
DBI action as explained in section 6.2. Since the fluctuation of decouples from every other 
fluctuation we can write down an effective Lagrangian for this fluctuation. This Lagrangian 
is given by 
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Using the results for the symmetrized traces ( ^.28 ), the effective Lagrangian simpUfies to 
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To simpHfy the equation of motion for the fluctuation (see ( 6.15| )) we define 
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